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THE MC-DAGUM DISTRIBUTION AND ITS STATISTICAL 
PROPERTIES WITH APPLICATIONS 


BRODERICK O. OLUYEDE, SASITH RAJASOORIYA 


ABSTRACT. In this paper, a new class of distributions called Mc-Dagum dis- 
tribution is proposed. This class of distributions contains several distributions 
such as beta-Dagum, beta-Burr III, beta-Fisk, Dagum, Burr III and Fisk dis- 
tributions as special cases. The hazard function, reverse hazard function, mo- 
ments, mean residual life function, Renyi entropy and Fisher information are 
obtained. Lorenz, Bonferroni and Zenga curves are derived. Maximum likeli- 
hood estimates of the model parameters and numerical examples are given to 
illustrate the usefulness of the proposed class of distributions. 


1. INTRODUCTION 


Camilo Dagum [1] proposed the Dagum distribution to fit empirical income 
and wealth data, and also accommodate heavy tailed models. Dagum distribution 
has both Type-I and Type-II specification, where Type-I is the three parameter 
specifications and Type-II deal with four parameter specification. The cumulative 
distribution function (cdf) and probability density function (pdf) of the three- 
parameter Dagum distribution are given by 


(1) G(z; ,9,8) = (1-- 078), 
and 
(2) g (1; 4,8,8) = 86x51 (1-- Ax?) ^ ^, for A, ô, 8 > 0, 


respectively, where À is a scale parameter, and ô and £8 are shape parameters. 
Dagum refers to his model as the generalized logistic-Burr distribution. Actually, 
when 8 = 1, Dagum distribution is referred to as the log-logistic distribution [13]. 
'The most popular Burr distributions are Burr-XLL distribution, often called Burr 
distribution with cdf, 


(3) F(zó,8)—1—(1-27) ^, forz 0,8, 8 » 0, 

and more importantly the Burr-III distribution with cdf 

(4) F(z;0,8) = (1-279) ^, for z » 0 and ô, 8 > 0. 
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Dagum distribution is a Burr III distribution with an additional scale parameter 
(A), [13]. The kt” moment of Dagum distribution is given by 


ô 
for ô > k, 4,6, B > 0, where B(.,.) is the beta function, (by setting t = (1+Aa~°)~!). 
The q'^ percentile of the Dagum distribution is 
(6) x(q) =A3(q? -1)*. 


Additional results on the Dagum distribution, including estimation of the parame- 
ters for censored data and properties of the beta-Dagum distribution can be seen 
in McDonald and Xu [14], Domma and Condino [2] and Domma et al. [3]. 


(5) E(X*) = BB (8+ 51-3). 


1.1. Basic Utility Notions. Some useful functions that are employed in subse- 
quent sections are given below. The gamma and digamma functions are given by 


Fee fe dandv(s)— LM respectively, where I” (v) = fj" t"-! (log t)e™ dt 
is the first derivative of the gamma function. The n'^-order derivative formula of 
gamma function is given by: I) (s) = if z571 (log z)” exp(—z) dz. The lower in- 


complete gamma function and the upper incomplete gamma function are 


(7) risa) = f t5 le dt and r(s,2) = f fae dt, 
0 x 





respectively. 


1.2. Mc-Donald Generalized Distribution. Consider an arbitrary parent cdf 
G(x). The pdf f(x) of the Mc-Donald generalized distribution is given by 


tah esc Ge IGT Gal, so poco Strands SO. 
B(a, 6) 

Note that g(x) is the pdf of parent distribution , g(a) = dG(x)/dx, and a,b and c 
are additional shape parameters. Introduction of this additional shape parameters 
is specially to introduce skewness. Also, this allows us to vary tail weight. Cordeiro 
et al. [6] presented results on the McNormal distribution. Marciano et al. [4] 
obtained the statistical properties of the Mc-I° distribution and applied the model 
to reliability data. It is important to note that for c—1 we obtain a sub-model 
of this generalization which is a beta-generalization and for a = 1, we have the 
Kumaraswamy (Kum) generalized distributions. If a random variable X has the 
pdf above, we write X ~ Mc-G(a,b,c). The cdf for this generalized distribution is 
given by 





1 G(x)* 
(9) F(z;a,b,c) = Ic. (a, b) = Sc wH = v)" dw, 
where L,(a,b) = B(a,b)~! [98V w^-1(1— w)^-1dw denotes the incomplete beta 
function ratio (Gradshteyn and Ryzhik, [8]). The same equation can be expressed 
as follows: 
G (z)^* 


(10) F(z;a,b,c) = aB (a,b) [2/1 (a, 1 — b; a + 1; G(x)?)], 
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where 





ET 1 po-1 1-— c—b—1 
(11) 2F; (a,b; c; x) = B (b, c — b) P (1- t) di 


0 (1 = tz)? 
is the well known hypergeometric function (Gradshteyn and Ryzhik, [8]), and 
B(a, b) = 2. 

An important motivation for the development of the McDonald-Dagum (Mc-D) 
distribution is the benefit of this class in its ability to fit skewed data that cannot 
properly be fitted in many other existing distributions. Mc-G family of densities 
allows for higher levels of flexibility of its tails and has a lot of applications in various 
fields. See Cordeiro el al. [5] for additional results on generalized distributions and 
McDonald Weibull distribution [7]. 

The hazard and reverse hazard functions for the Mc-G distribution are given by 


f(wja,b,c) ww) (x) {1 - e 











(12) hp (x) = "Em F(x: a, b, c) = B (a, b) fi = Tage (a, b)} ) 
and 

í ac—1 T — G(x b—1 
T we (a) = ith. cg (2) GA? (x) (1 - GG) 


F (a5 a, b, c) B (a, b) Iac(z) (a, b) 


respectively, for a > 0, b > 0, and c > 0. 

The outline of this paper is as follows. In section 2, the Mc-Dagum distribution 
and related family of distributions are introduced. The series expansion for the 
density, hazard and reverse hazard functions, and other properties are presented in 
section 3. Section 4 contains the moments, Lorenz, Bonferroni and Zenga curves. 
Section 5 deals with measures of uncertainty, including Renyi entropy. Section 6 
contains maximum likelihood estimates of the model parameters. Fisher informa- 
tion and asymptotic confidence intervals are presented in section 7, followed by 
applications in section 8. Some concluding remarks are given in section 9. 


2. Mc-DAGUM DISTRIBUTION 


In this section, the new class of distributions, called McDonald-Dagum (Mc-D) 
distribution is introduced. Considering the properties and some useful features 
of both Dagum and Mc-Donald distributions, a broad range of generalization is 
possible by combining these distributions. The new class of distributions possess 
capabilities widely applicable in several areas including economics, finance and re- 
liability. 

Now, combining the Mc-G and Dagum distributions, we obtain the pdf of the 
Mc-Dagum distribution as follows: 


R - cBAóx- 5-1 —§\—Bac-1 —§\—cB b-1 
f(x; A, 6, B, a, b, c) a ar AC ae eh ) j- (1+ A2 ) | $ 
for a,b,c, 4, 8,6 > 0. The cdf of this new distribution is given by 
(14) F(x) = Lapas). (a, b), 


where 
I, (a,b) = BED To: q-1(1— w)*-idu 


4 RODERICK O. OLUYEDE, SASITH RAJASOORIYA 


is the incomplete beta function. The cdf can also be written as follows: 
(1+ ar78) "** 





Hr AE [2F; (a,1 — b; a + 1; (1 + A 79)799)], 
where 

b—1 c—b—1 
(16) 2F; (a,b; c; x) BO 5 J y i —— dy, 


is the well-known hypergeometric function, (Gradshteyn and Ryzhik,[8]). 
The hazard and reverse hazard functions are given by 


c 75-1 qz -Bac-lp, — 28) 688-2} 
hp (z;>,6,B,a,b,c) = 2% (1&327*) [iiras 


B(a,b) E agen] 





(1+às78) 
and 


cBAx- 9-71 1+Ax75 ni um 1—-(14-Ax-? ce 
TF (x; A, ô, B, a, b, c) 2 B( | a,b) ) 





aba 4 9-8) - Bel 
for a > 0,b > 0,c > 0,à > 0,8 > 0,6 > 0, respectively. 
3. EXPANSION OF DISTRIBUTION 


Note that for |w| < 1, and b > 0, (1 — w)®™! = 2558 CRW i, Therefore, the 
cdf can be expanded to obtain 





po pO) 2 (—1)/T (b 
Paassen = aly o vy eas 
(17) m3: 0 5005) 
j=0 
where p; = BER eee Similarly, the pdf is given by 
(18) djs Cero ey 
j=0 


If b is an integer, then the summation in equations (17) and (18) stops at b — 1. 
If c = 1, we have a finite mixture of Dagum distribution with pdf 
b—1 
f(x; A, 8,0, a, b) Tz X pj 9(2; B (a + j), À, ô). 
j=0 
The graphs below are the pdf, cdf and hazard function of the Mc-Dagum distri- 
bution for different values of parameters A, ô, 8,a,b, and c. The graphs of hazard 
function shows various shapes including decreasing, unimodal shapes for the select- 
ed values of the model parameters. 
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3.1. Submodels. With this generalization, we have several sub-models that can 
be obtained with specific values of parameters A, 8, ô, a,b and c. 


(1) When c = 1, the Mc-Dagum distribution is the beta-Dagum Distribution, 
with the density given by: 


BAóz- 5-1 
B(a, b) 
for x >0,A>0, 8 > 0, ô > 0, a > 0, and b > 0. See Domma and Condino 


[3] for additional details. 
(2) If a = b = c = 1, we have the Dagum distribution with the pdf, 


(20) fp (2; À, ô, B) = 86x71 (1 + Ar?) ^, 
for A, ô, B > 0. 
(3) If b = c= 1 and a > 0, then we have the Dagum distribution with param- 
eters Ba, and 6. The pdf is 


b—1 


(19  f(a;, B,6,a,b) = Geigy n d ee, 


—Ba-1 


(21) f (2; Ba, A, 8,) = Barden} (14 A73) T, 
for A, 6, B > 0. 
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(4) If a = c = 1 and b > 0,we have another Beta-Dagum distribution with 
parameters b, 8, ^, Ó and the pdf is given by 


(22) fap (2; A, ô, B, b) = 88392 (1-277) ^ [1 - (13275) ^] at 


for A, ô, B > 0. 
(5) Ifa = c = à = 1, then we have the beta-Burr III distribution with param- 
eters b, 3,6 and the pdf is given by 


(23) fap (1:0, 8,0) = bBóx -?-! (1 qug ors [ _ (14 iy p 


for 6,6, 8 > 0. 
(6) If c = 8 = 1, then we have the beta-Fisk distribution with parameters 
à, b, A,6 and the pdf is given by 


Ax?! 


(24) fer (x; à, ô, a,b) = “B (a,b) 


(i632) fi- 0427) 7], 
for a,b, 4,0 > 0. 


3.2. Kum-Dagum Distribution. Kumaraswamy [11] in his paper (1980) pro- 
posed a two-parameter distribution (Kumaraswamy distribution) defined in (0, 1). 
Here we will refer to it as Kum distribution. Its cdf is given by: 


F(x;a;b) 21—(1—a?)') 2x€(0,1),a>0,b>0. 


The parameters a and b are the shape parameters. The Kum distribution has the 
pdf given by: 


f(z;a,b) aba" !(1—27)!, 2 €(0,1),a>0,b>0. 


Note that the Kumaraswamy distribution can be derived from the beta distribu- 
tion. Combining cdf of Kumaraswamy distribution with the Dagum distribution 
discussed in earlier, we obtain Kum-Dagum distribution with the cdf and pdf given 
by 

— Ba? 
(r)-1-— l- (1+ Ax?) $ | 


Kum 


and 


—Ba- —BaN —£— 
ese) = ipit oe] (1 = [raat 


for a,b, 8,A,6 > 0, respectively. We do not study the properties of the Kum-Dagum 
distribution in this paper. 


4. MOMENTS AND INCOME INEQUALITY MEASURES 


In this section, we present moments, Lorenz, Bonferroni and Zenga curves for 
the Mc-Dagum distribution. Income distribution and its variation is an important 
concern for economists. We use the results presented earlier, which was obtained 
by expanding the pdf. 
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4.1. Moments. We can derive the k*” moment of a Mc-D distribution using prop- 
erties of the mixture distribution. The k*” raw or non-central moments are given 
by 

(25) 


= b—1 
E(X*) - fox gh Bde = a IE te (1 = (1 + A07#)"*) di 
b—1 


" for ak FA (1985) a (23-5) 7) ds. 
Now let y~! = (1+ Aw~*) then z = (1— y) * (Ay)*, and we have 
E(X*) = ggf fo (1 — y) E Oy) tyt — ye) Ty 


Using the fact that (1—y*)^7! = 577. n (P), and for pj = ae. 


we obtain 


Be oo (-1)/ (b ups CES 
EQ) = de Xo PEIN Jo vi etti — y)dy 


= eh Ds FERRO ee) el - b 


= Ye, ACE Bola 4. j) + 8,1— 5), 








for ô > k. 
The mean residual life function (MRLF) denoted by u(x; 2, 8, ô, a, b, c) = u(x) is 
given by 
(26) 
we) = EJX —z|X > a] 
— E(X)- E(X|X Ev) 
5e 1— F(x) 
me p; Beas) 
j=0 5 





B(Be(la+j)+#,1-8)- 55 pj step B((14 ArT?) l; Be(a+j)+8,1— 
1— ee SC ASUEI à) 





ole 
— 








4.2. Lorenz, Bonferroni and Zenga Curves. Lorenz, Bonferroni and Zenga 
curves are the most widely used inequality measures in income and wealth distri- 
bution, [9]. In this section, we will derive Lorenz, Bonferroni and Zenga curves for 
the Mc-Dagum distribution. The Lorenz, Bonferroni and Zenga curves are defined 
by 








[rtf (at [i tf (Oat 

L muU IET d B(F o AONE 
and A(x) = 1-419 , respectively, where u^ (x) = Jo tf Fe age and ut(x) = fe fee = 
E reii are the lower and upper means. For Mc-Dagum distribution, using 


these results, we obtain the curves. Lorenz and Bonferroni curves for Mc-Dagum 
distribution are given by 


Moo PiBc(a + j)A3B B((1 + Axr?) t; Bela + j) + i 1- i) 


TOJ = 
VID). EnS) Eo Di Bela + j)A* B(Be(a + j) + $,1 - $) 





? 


and 
(28) 


- ;Dc(a NAS z79)-1: Bela . +4 201 
B(Fa(z;0)) = jo PjBc(a + 3)^* B((1 + Az?) =t; Bela 3) + $,1— 3) 





SEGE A Bola 2,0 Se orba A B(Bea + 1 D 
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respectively, where O = (A, 8,6, a, b, c). Zenga curve for the Mc-Dagum distribution 
is given by 





: - Q- F(z))E[X|X xz] 
(29) A(r;,Qg) =1 F(z)|E(X)-E(X|X €z)] 


where E[X|X € a] 2 Xo p; —— p —— folat i)a? * B((1+Ax~8)—1; Bc(a-- 5) - 1,1— 1), E(X) = 
7? o pj Bc(a--j n ; oo - 
Mosi?" B(Bc(a + j) + 4,1— 3), and F(a) = IZo pG (2; A, Bela + j), 6). 


5. MEASURES OF UNCERTAINTY 


_ In this section, we discuss the Renyi entropy [15], Shannon entropy [16] and 
B-entropy for the Mc-Dagum distribution. 


5.1. Renyi and Shannon Entropy. For a pdf f(x), the Renyi entropy is given 
by 








l oo 

(30) Ar(f)= : ET e Pius). s 7» 0, and s Z 1. 
= 0 

As s — 1, we obtain the Shanon entropy. Note that, 

a E (cBAS) x75275 _§\—Bacs—s —§\—cB bs—s 
P (cba) [9 ó— —ô\— — —5\—cBy}bs— 
f(x)dv = yp 99 S(1 + Ax PRESS cpu Aue cPjbs-s dy 
] re = Eos], (haart) eem (1 ae) 








cpad)’ ; sd—s, 78S 4 Bacsts Ei c\sb—s s— 


Using the fact that, for |w| < 1, (1 —w)^71 — 3g CELO, and setting y — 
(l-FAz-9)-1, so that z7? = y = PU: and Aóx-?-!dz = y-?dy, we obtain 

(31) 

D Idee ACH CE rBaeske- 2-1 (4 — yBejsb-e(1 — yet $-$-1dy 











= Bs(a,b) 
——. (e8A8)* Att g-sb-s 84 Eo 1 4s8—1_] (— DA (sb— s--1)y^ €? 
- BG) oy Bacs— $ (1- y)? 85$ xu 0 — Fb sHI-j)jl dy 











_ (epas) atta- sre, EVM sb=s+1) p(gc; 


S 
= Bs (a,b) j=0 T(sb—s+1—))9! + bacs — $ + 1,5 $ — D). 


Therefore, Renyi entropy for the Mc-Dagum distribution is 





log [ore $ a —1))I'(sb—s+1) 


(32) Ar(f) TEF B*(a,b) FG SIS 


j=0 


. s 1 s 1 
B (5s + Gacs 57 P? "5 5) | 


for s > 0 and s Z 1. If bs — s is a positive integer, then the sum in the Renyi 
entropy stops at bs — s. 
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5.2. B-entropy. We also obtain B-entropy for the Mc-Dagum distribution. Note 
that B-entropy is given [o 


Hie -f f? (a (ode, if B £ 1, and B > 0. 


pe 
If 6 — 1, then Hg(f) = E[-log(f(X))] is the Shannon entropy. 
For B #1, and B 0, B-entropy for the Mc-Dagum distribution is 
1 f (cBAS)PAL +3 -3-A 5S (-1) T(Bb — 8 +1) 
po BP (a, b) Spt (b= Bop tit 











X B (sj + Bach A J| 





6. MAXIMUM LIKELIHOOD ESTIMATES 


Let O = (A, 8,6, a, b, c)T. In order to estimate the parameters A, 8, ô, a,b and c of 
the Mc-Dagum distribution, we use the method of maximum likelihood estimation. 
Let £1, Xo, ....... ,Zn be a random sample from f(x; A, 8,6, a, b, c). The log-likelihood 
function L(A, 8,0, a, b, c) is: 


L(A,B8,0,a,b,c) =  mlog(c) + nlog(B) + nlog(À) + nlog(0) — nlogB(a, b) 


— (6 4- 1) V 7 logz; — (Bac +1) X` log[1 + Az; ?] 


i=l i=l 


n 
- (b — 1) V ^ log[1 — (1 + A27?) °]. 
i=1 
Differentiating L(A, 8,6,a,b,c) with respect to each parameter A, 8,ô,a,b and c 
and setting the result equals to zero, we obtain maximum likelihood estimates. 
'The partial derivatives of L with respect to each parameter or the score function is 
given by: 














OL OL OL OL OL OL 
(33) v.) = (zai ar aE or ac) 
where 
hos 
n —ó n a n —d)—cB-1,,—-6 
cB(1 + Arz’) 8-12; 
Bac > 2s te-05, ; 
m —d)\_¢ 7 
-se ( J eum Um c Palade ra 
OL n ™ (14+ Az; )-*Plog(1 + Av; ?) 
35) — log(1+ À +c(b-1 : a 
( ) 9B B -aX og( Xu 2) c( 2. [1 o (1 "E Azz?) <$] 
OL a; ^elog(z;) 
> = l i) +A +1 n Ee 
86 3 - tow) Bae Mu =) 
X. MÀ ys ley) 
36 — xX , 
(36) cB(b Dx Ger] 
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(37) ac = -n(éa) = Wla +0) = Be log(t + daz", 
OL us 


(38) = —n[v(b) — yla + b)] + X` log[1 — (1 + Az; ?) 2], 


i=l 


where w/(.) is digamma function (x) = “logI'(x) = LE , and 


OL n 3 -ó = "(1+ A2; ?)-9810g(1 + Ax; ?) 
(39) 57-7 Bao Mi + A2; ?) + (b pe Hen us 








The MLE of the parameters A, 8, ô, a, b and c, say d B, ô, a, b and ĉ are obtained 
by solving the following equations, oL on ob aL oL 2L 0. There is no 











closed form solution to these equations, so numerical technique such as Newton 
Raphson method must be applied. 


7. FISHERS INFORMATION MATRIX 


To obtain the Fishers information matrix (FIM), we derive the second partial 
derivatives and cross partial derivatives with respect to each parameter A, 8, ô, a, 
b, and c as follows: 

OL -n i a; (6-1) Y eBay? A; PPAP — eB — 1] 
[1 — A; ^p 





, 








3 


^L n m n of _ep—1 [eBlog(A;) — 1 + A; ^^] 
(41) - 2? -+ (b-1) 5o; ?cA; ae 





OL 2 2; ^log(zi 2 ae 
(42) = (— fac nyo og (i) + (b 1) X` cfiz; A; 4 1og(x;)Bi, 
i=1 i i=1 
1—Az; cBAT; 5 —Az; A; 4A; ^ cBA; 0 1 - (14-3; 2) 78 

















where B; — i ap 
(43) DL —Be > za 

i —5\)—cB-1,,-6 
(44) ns = i an 2 : = : 
Wb gaT AE (ü- np A Ee LR AUT, 
(46) a -x Mel a+ n à PER 


ƏL z Z Y xill — cBlogA; — Ar? 
I uya- 9; g 


i [are 





’ 
i=l 
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Ax ôl ; 
where C; = 25729857 IC 























Ai ? 
OL i j 
4 =e ` log(1 + ee 
(48) dias 7-2; toatl +279) 
(49) PL 5 e(1 + Ax; 9)- 8 log(1 + Ax; ?) 
080b — [L— (1222; ))-98] | ' 
"A, logA;|cBlogA; + A; — 1] 
50 —a logA; + (b— 1 t L ; 
E AD, Oe 
27, "E n 
(51) T = > + A(Bac + 1) si) (logi A; ^! D 
i=1 
ais Acx; 5(14-2; 9)7 l—Ax; Hee aes B+ ar 5(14-32; 2)- 99-1] 
so, 
(52) OL A8 a a, loga; 
—— MeN i 
060a — (1+ Az, ?) 
OL az; 9(1 + Ax; ?)-99-1ogz; 


M o ' 





o L 2 7 A; 71, -logz;[cBlogA; + A; ? — 1] 
54) == =À F; — AB(b—1 4 t ; 
(54) c 2? BC ) th A; 9] 


lo Z4 
where F; = 2 109% 


























(1+Ax;°) 
65) Shan le eap - TD - quay + FO], 
(56) A = jou ro- E, 
(57) = = A2, log(1 + àx7°), 
Gm =n lee t- TES - wort Te. 
OL XB(-c-2z?)- *log(1-- Az; ?) 
em Obüc - x L-art] — ' 
and 
2 = n a7 —cB o aT) 2 


[Ies (poa ^) «ep 
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Fisher information matrix (FIM) for the Mc-D distribution is: 


Ixx Tye dxs Da Do Ixe 
lg. gg Ips Iga Ige Ipc 
| ls. deg los sa Isp sc 
Tes " xj 4 9 B Tar lag Las Loa Lav Lac i 
Ioa Lop Tas Toa To lw 
Ley Leg Les Lea Io Lee 


where ha = —E OL neat Je = -E = . The elements of the 6 X 6 matrix 


I(A, 6,6, a,b,c) can be approximated by the elements of the observed information 
matrix, where 





= In(®). 


(ie r| A -PL 


50,00; | ~“ 00,00; 


7.1. Asymptotic Confidence Intervals. In this section, we present the asymp- 
totic confidence intervals for the parameters of the Mc-Dagum distribution. The 
expectations in the Fisher Information Matrix (FIM) can be obtained numerically. 
Applying the usual large sample approximation, MLE of O, that is, Ó is approxi- 
mately Ne(O, I, 1 (O)), where I, (O) is the 6X6 observed information matrix. Under 
the condition that the parameters are in the interior of the parameter space but 
not on the boundary, the asymptotic distribution of /n(Ó — O) is N«(O, I-!(O)), 
where 7(0) = Jim. n~! (O) and 


Da Te hè Da Do Le 

Ig, gg Iss Iga gp gc 

= Is, Isg Iss Isa Isp Lee 
I,(8) RM Tar lag Las Loa lab Iac 
Isa Lop Tes Tba To Toe 

Ley Leg Les Lea Ie» Loe 


The multivariate normal distribution with mean vector (0,0, 0,0,0,0)7 and co- 
variance matrix /,(O) can be used to construct confidence intervals for the model 
parameters. That is, the approximate 100(1 — 7)% two-sided confidence intervals 
for A, B, 6, a, b and c are given by: 


At Za 1: (6), B+ Za I7}(6), St Za 15; (8), a+ Za Iz (0), 


b+ Zay Ip (8) and @+ Zy Iz (6), 


respectively, where Z a is the upper a percentile of a standard normal distribution. 

We can use the likelihood ratio (LR) test to compare the fit of the Mc-D distri- 
bution with its sub-models for a given data set. Note that, to test b = c = 1, the 
LR statistic is 


wW = 2[In(L(B, ô, Pr à, b, é)) m In(L(B, ô, À, a, 1, 1))], 




















where B, ó, Â, à, b, and 6, are the unrestricted estimates, and B, ô, À, and à are the 
restricted estimates. The LR test rejects the null hypothesis if w > XS. where x2 
denote the upper 100d% point of the x? distribution with 2 degrees of freedom. 
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TABLE 1. Estimation of models for failure times data. 









































Estimates Statistics 
Distribution B ô À a b c —2In(L) AIC AICC BIC 
Mc-D 4.1652 2.7600 106505 0.01904 45431 4.1652 26.0 38.0 148.5 142.3 
(5.3470) — (0.3228)  (2.966E-6) (0.01999) (0.04933) (5.3470) 
B-D 7.2221 8.8376 3.5995 0.6595 0.04258 - 47.1 57.1 163.8 160.6 
(2.9051) (0.5793) (6.9162) ^ (0.8025) (0.01171) : 
B-BHI 0.8308 0.09752 - 290.58 172.49 - 31.3 39.3 143.3 142.1 
(0.03511) (0.01684) - (0.2349) ^ (0.1349) - 
D 1.1568 1.5295 77.9568 - - - 31.6 37.6 139.8 139.7 
(0.6532) — (0.3450) (123.96) - = 
BHI 1.3821 1.2237 - - - 36.0 40.0 141.0 141.5 
(0.04983) (0.03518) - ^ - z 
Mc-W 8.8539 0.01463 0.006759 362.23 13.2979 26.0 36.0 142.5 139.5 
(0.4082) (0.001952) (0.009667) (12.1191) (18.8857) 





8. NUMERICAL APPLICATION: Mc-DAGUM AND SUB-DISTRIBUTIONS 


In this section, applications based on real data, as well as comparison of the Mc- 
Dagum distribution with its sub-models are presented. We provide examples to 
illustrate the flexibility of the Mc-Dagum distribution in contrast to other models 
for data modeling. These data sets are modeled by the Mc-D distribution and com- 
pared with the corresponding sub-models, including the Dagum distribution. The 
first data set (Lawless [12]) represent the failure times, in minutes, of 15 electronic 
components in an accelerated life test and they are as follows: 1.4, 5.1, 6.3, 10.8, 
12.1, 18.5, 19.7, 22.2, 23, 30.6, 37.3, 46.3, 53.9, 59.8, 66.2. The second data set 
presented in Tables 2 and 3 represents the salaries (in dollars) of 818 professional 
baseball players for the year 2008. We fit Mc-Dagum(Mc-D), beta-Dagum (B-D), 
beta-Burr III (B-BIII), Dagum (D) and McDonald Weibull (Mc-W) distributions 
to these data using the method of maximum likelihood estimation. The MLEs of 
the parameters (with standard errors in parenthesis), Akaike Information Criterion 
(AIC), Bayesian Information Criterion (BIC) for the fitted models are presented in 
Tables 1 and 4. The Mc-Weibull pdf [7] is given by 
cô Ax?! — (Ax)? 

B(a,b) 
for A, 6,a,b,c > 0 and x > 0. 

The maximum likelihood estimates (MLEs) of the parameters are computed 
by maximizing the objective function via the subroutine NLMIXED in SAS. The 
estimated values of the parameters (standard error in parenthesis), -2log-likelihood 
statistic, Akaike Information Criterion, AIC = 2k — 21n(L), Bayesian Information 
Criterion, BIC = kln(n) — 21n(L), and Consistent Akaike Information Criterion, 
AICC = AIC + 2***) where L = L(Ó) is the value of the likelihood function 
evaluated at the parameter estimates, n is the number of observations, and k is 
the number of estimated parameters are presented in Tables 1 and 4 for the Mc-D 
distribution and its sub-distributions. 

For the failure time data set, the likelihood ratio statistics for the test of the 
hypotheses Ho : B-D against Ha : Mc-D and Ho : BIII against Ha : Mc— D are 
21.1 (p-value« 0.00001) and 10 (p-value « 0.05), respectively. Also, the likelihood 
ratio statistic for the test of Ho : B-BIII against Ha : Mc-D is 5.3 (p-value=0.07065). 
Consequently, we reject the null hypothesis in favor of the Mc-D distribution and 
conclude that the Mc-D distributions are significantly better than the Dagum and 
Beta-Dagum distributions based on the LR statistic. The values of the statistics 


(61) f(z;A,,a,b, c) = (Seer eh ae Deed 


, 
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TABLE 2. Baseball Player Salary (x109) 
0.403 1.75 6 0.4345 4.956237 5.5 0.78 0.475 L5 
11.660000 13 0.4186 2.5 0.406 13.1 1.775 2.7875 | 04 
0.8 6.25 0.404 1 1.0825 0.42245 9 4.05 0.4 
8.75 1.75 5.9 1.75 4 0.4551 3.125 0.975 5.5 
1.5 5 1.5 1.7 10 15 0.4073 14 8 
6.25 0.441 3.65 2 0.800002 33 0.4 1.98125 — 0.424 
0.5 1.5 0.4363 3.5 1 15.285714 1.25 3.666666 0.75 
0.401 5.5 0.4142 0.4275 0.403 5.4 0.4115 7 0.4 
7.5 0.4 1.95 19 0.4 20.625 0.5 0.675 0.452 
3.05 5 4.766666 5.5 7 0.432975 — 0.4044 8.25 0.445 
3.5 4 0.4275 0.75 0.414 5 0.4324 — 0.4333 2.8 
0.425 6.25 10 2.8 6.925 0.4 0.4300 1.255 0.475 
0.7125 7 10 0.75 4.65 0.4 0.4 0.4337 — 0.425 
0.43 11 7.05 14 3.25 0.405 0.41 0.75 0.5 
0.4115 9 0.415 0.439 6 0.41 11 3.25 2.95 
2.535 0.43 1.625 0.61 0.95 0.41 0.675 ^ 0.4104 5 
2.525 12.5 1.055 1.5 5 4 1.6875 1.000018 — 0.4 
0.4225 0.8 9.375 2.6 4.75 0.4 0.4 6 3.325 
0.4 4.2 0.75 0.449 1.6625 0.42 0.4005 0.55 
1.45 0.4215 2 0.5 0.4 11.5 4.625 24 
3.6 7 2 14.811414 0.4 0.4 3.5 1.8 
0.4225 1.625 0.401 8 0.426 5.8 1 1.29 
2 0.475 3.75 0.425 4.25 2.8 0.404 8 
5 0.4075 0.8 0.4245 0.415 0.41 2.9 4.445 
0.412 3 5.5 1.2 0.4075 0.41 0.40175 — 0.4028 
6.5 13.25 16.6 1.4 12.75 3.9 2.75 0.4139 
1.95 13 0.415 2.9 0.447 5 0.401 1.5 
2.5 0.575 0.75 0.421 2.45 0.41 8 1.35 
0.4375 0.45 9.875 0.4155 1.3 0.41 0.44 0.4301 
0.4325 2.2 0.402 0.411 0.75 1 0.65 0.435 
0.42 0.4 0.401 6.083333 0.44 4 4.75 3.2 
0.425 16.65 0.408 2.25 0.45 13.5 0.4 0.4159 
1.2375 17 0.75 1.3 0.422 0.4 0.41 0.4 
0.425 0.575 0.403 3.6 0.4125 0.4 0.475 3.75 
3.375 0.5 0.407 4.25 0.4275 3.3 6.5 0.75 
1 3.5 0.415 0.457 1.7 0.41 8 6.2 
3.45 18.75 0.404 3.75 2.5 0.41 0.4017 8.06 
15.5 0.44 4.5 12 7.666666 1.1 9.6 0.40807 
1.85 1:3 3.75 3.275 0.525 0.405 0.455 0.40662 
2.825 14 1 4 0.4 5.475 — 0.40175 0.41 
10 0.44 0.4 114 0.4375 0.8 0.40125 — 0.43468 
0.45 1 12.5 2.7 1.475 0.5 0.4 1.615 
8.333666 10 0.4 0.425 0.4661 0.65 1 1 
2.885 0.52 14.883049 0.405 1.15 0.46 0.401 0.555 
15 6 1 1.8 2.75 1.15 0.404 0.4055 
8.7 11.5 0.435 1 0.43 1.635 0.48 0.401 
0.405 0.41 3.5 3.575 10.5 2 9.6 0.8 
0.41 0.75 10 0.437 0.4 5 18.5 0.5 
0.8 0.401 6.3 0.64 11.6 4.35 0.419 3.2 
0.4 5.6 2.2 5 11.25 0.405 0.405 0.40473 





AIC, AICC and BIC shows that the Mc-D distribution gives smaller values than 
B-D, B-BIII, D, BIII, and compares favorably with the Mc-Weibull distribution. 

For the baseball player salary data set, the likelihood ratio statistics for the test 
of the hypotheses Hp : D against Ha : Mc-D and Ho : B-D against Ha : Mc-D are 
519.4 (p-value« 0.00001) and 14.6 (p-value « 0.0002), respectively. Consequently, 
we reject the null hypothesis in favor of the Mc-D distribution and conclude that 
the Mc-D distribution is significantly better than the Dagum and beta-Dagum 


distributions based on the LR statistic. A closer look at the values of the statistics 


AIC, AICC and BIC shows that the Mc-D distribution gives smaller values than 
B-D, B-BIII, D, BIII, and Mc-Weibull distributions. 
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TABLE 3. (Continued) Baseball Player Salary (x 10°) 
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0.415 12 2.8 0.4 0.44 15 9.5 0.65 
1.4 0.4025 0.403 0.41 0.43 7.166666 7.75 12.868892 
0.4 4.5 4.25 0.46 4 12 13.4 0.41482 
0.4 2.65 0.4375 0.403 0.95 2.333333 2.3 0.41176 
6.35 0.45 18.971596 10 0.44 6.5 1.4 12 
11.5 0.4 0.4 5.775 0.435 12.083333 0.8225 0.41631 
0.41 6.25 4.6 8.5 1.152 2.5 2.05 0.41089 
0.4 0.55 2.095 15 1.15 8.5 0.413 0.401 
7.166666 1.225 0.75 18 19.243682 4.25 0.405 1 
0.405 0.405 7 1.6 6.25 0.475 2 1.85 
4 13 2.7 0.401 2.625 1.625 0.455 13.054526 
0.65 1.325 0.825 0.465 2.8 0.835 3.8 2.5 
1.5 0.55 1.3 10 2 11.285714 0.4 0.6 
0.445 10.125 2.275 0.405 12 3.125 0.405 2.4 
8 0.4175 0.75 0.45 1.6125 2.5 7.666666 0.75 
2.4 0.42 3.675 10.4 0.471 0.4 0.4 0.4087 
0.42 12.125 10 L3 2.6 0.4085 0.4 3.75 
0.435 0.418 0.735 0.435 0.401 2.425 2.25 1.45 
3.35 0.4 1.3 2 0.406 0.402 0.418 14.25 
0.8 2.425 0.4 3.665 0.575 0.4015 0.425 2.59 
8 5.375 0.4 1.1 0.4 2.15 0.42 0.4139 
0.41 2 0.5 3.25 2.2375 2.2 12.25 0.64 
0.85 0.4 3.5 3.8 12 0.4015 18 0.4118 
8 14 0.4 0.475 6 2.3 0.5 0.4144 
0.41 8.5 1:9. 5 2.25 0.825 0.4 0.4052 
2.4 0.4 0.42 4.75 0.925 1.5 9.85 1.9 
2.75 0.4 0.95 0.465 6.125 0.4 2.825 0.4194 
0.4 1.5 2.25 9.166666 0.4015 0.4 0.85 
5 0.418 1.4 0.475 18.876139 7.05 13.302583 7.95 
4.5 0.42 0.66 5 4.9 0.4135 0.825 0.4037 
2.5 5.1875 0.41 1.825 14 2.5 2.5 0.4023 
1.5 2 0.4 0.405 0.4095 2.5 0.4 6.4 
0.5 2.25 0.4 3.1 Lz. 0.75 0.4 11.625 
7.8 0.44 2.4625 7.5 10.5 0.4115 12.137 0.4 
11.166666 0.4375 2.4 3.364877 7.75 0.4145 6.5 12 
0.4 3.5 0.4125 0.467 0.403075 0.4135 0.405 1.1 
0.476 0.4161 0.44 0.404 1.25 6.25 0.95 0.4014 
14 3.35 0.4 0.437 16.5 3.2 7.4375 1.015 
0.4495 0.4167 0.404 12.433333 1.4 1.875 3.7 4.6875 
0.4 2.9375 0.75 0.465 6 0.4 0.411 1.6 
0.55 2.5 5.5 1.25 0.432575 7.4 0.5 2 
1.5 0.4203 2.225 3.9 0.4033 1.3 3.3125 1.9 
0.4 0.4 0.4 0.40075 13 0.4148 0.4 2.6 
1 5.5 5.35 2.35 0.414 0.4299 7.5 0.431 
1.35 4 0.4 10 21.6 1.255 14.427326 0.4155 
12.5 0.4214 0.4 23.854494 3.75 0.75 0.4 8 
0.414 0.4461 1.55 15.217401 13 0.4 1 0.4155 
9.25 11.5 14.5 0.401 2.125 0.75 0.403 0.4 
0.4155 0.4038 1.25 0.4 6.55 0.4 0.43 8 
8.333333 3 0.75 0.4025 0.4 0.85 0.65 0.42 
TABLE 4. Estimation of models for baseball player salary data. 
Estimates Statistics 
Distribution B ò À a b € —2In(L) AIC AICC BIC 
Mc-D 47.9710 8.0983 0.000220 0.1053 0.06001 0.8718 2706.2 2718.2 2718.3 2746.4 
(9.9363) (0.05615) (0.000014) (0.01043) (0.002683) (0.1580) 
B-D 24.7269 8.2030 0.000317 0.1196 0.07237 - 2720.8 2730.8 2730.8 2754.3 
(1.0055) (0.02546) (0.000041) (0.01180) (0.003253) E 
B-BIII 8.4648 1.2126 - 0.1646 1.0218 - 3367.3 3375.8 3375.3 3394.1 
(2.4174) (0.06022) : (0.06231) — (0.09221) d 
D 70.078 1.0301 0.01163 - - - 3225.6 3231.6 3231.6 3245.7 
(34.9306) (0.03011) (0.005860) - - - 
BII 1.3821 1.2237 - - - - 3367.3 3371.3 3371.3 3380.7 
(0.04983) (0.03518) f > : 
Mc-W - 0.2562 68.7906 0.01455 4.9316 1315.50 3300.8 3310.8 3310.8 3334.3 
(0.01189) (26.4644) (0.003375) (0.1424) (354.12) 
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9. CONCLUDING REMARKS 


We have proposed and presented results on the mathematical and statistical 
properties of the Mc-Dagum distribution. This class of distribution contains a 
fairly large number of distributions with potential applications to a wide area of 
probability and statistics including income and lifetime data analysis. Properties of 
the class of Mc-Dagum distributions including the pdf, cdf, moment, hazard func- 
tion, reverse hazard function, inequality measures including Lorenz, Bonferroni and 
Zenga curves, Fisher information, Renyi entropy and (-entropy are derived. Esti- 
mation of the parameters of the models and applications to illustrate the usefulness 
of the distribution are also presented. 
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